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Abstract 



K.^ An external magnetic field is found to have strong effects on the electronic 

j~^ ' structure of carbon nanotubes. A field- induced metal-insulator transition is 

o' 

C^^ ' predicted for all pure nanotubes. In a weak field, nanotubes exhibit both large 



diamagnetic and paramagnetic responses which depend on the field direction, 



C^ ' the position of the Fermi energy, the helicity, and the size of the nanotube. 

a" 

Universal scalings are found for the susceptibility as functions of the Fermi 



energy, the temperature, and the size of nanotubes. These results are in 
agreement with experiments. 



'i_j ■ PACS numbers: 61.46+w, 36.40-hd 



Typeset using REVTgX 



The exciting discovery of carbon nanotubes [0] stimulated a large number of theoretical 
studies on their electronic properties. Both tight-binding 0J^ and first principles calcula- 
tions [^] predicted that nanotubes can be either metallic or semiconducting depending on 
their helicity and size. Several recent experiments demonstrated some unusual properties. 
Giant magnetoresistance and indications of a field-induced metal-insulator transition are 
found in transport measurements |^. Large diamagnetic susceptibilities are found for a 
magnetic field both perpendicular and parallel to the tube axis PJ^. 

The magnetic properties of nanotubes were studied by Ajiki and Ando |^ using the 
k • p perturbation method However, they found that the susceptibility xi. (fo^' the field H 
perpendicular to the tube axis z) is three orders of magnitude larger than x\\ (when H || z). 
This finding disagrees with experiments, where it is found that the two are comparable 0. 
The k ■ p calculation is valid only if the Fermi energy is at the center of the band (half 
filling), and it only provides information about the band structure near the Fermi energy. 
The orbital magnetism depends on the total band energy which requires the calculation of 
the total band structure. In addition, in real materials it is likely that the system is not 
exactly at the half filling (corresponds to the case of a finite carrier density in graphite). 
Thus, a calculation which includes the field dependence of the complete tt band is necessary 
to understand the magnetic response of the nanotube. 

In this letter we report the results of such a calculation using the tight-binding model and 
the London approximation. We obtained the following results. 1) A magnetic field induced 
metal-insulator transition is predicted for all pure nanotubes, the transition depends on 
the helicity, the radius i?, and the magnetic field direction. 2) The weak-field magnetic 
susceptibility is large and increases linearly with the size of nanotubes, x ~ -R; it can be 
either diamagnetic or paramagnetic and is sensitive to the position of the Fermi energy ep- 
3) Associate with each nanotube is a unique energy scale Aq; the scaled susceptibihty xl R is 
found to be a universal function of e^/Ao and ksT /Aq for each family of nanotubes. 4) For 
typical nanotubes, x ~ — 250xl0^^cgs/mole, and |x| decreases with increasing T. These 
results are in agreement with recent experiments 



We use the nearest-neighbor tight-binding hamiltonian to calculate the band structure 
formed by the vr orbitals. This hamiltonian has been shown to be an excellent approxima- 
tion for calculating the electronic structure of fuUerene-related materials such as the large 
fuUerene molecules P, the solid fuUerite and fuUerides [0, and the nanotubes 0. Including 
the effect of a magnetic field in such a model is straightforward in the London approximation, 
which has been used successfully for studying the ring current and the magnetic response 
of Ceo and C70 molecules |jll|. The symmetry of nanotubes have bees studied by several 
groups 0. We follow the elegant approach of White et al. |^. The structure of a nanotube 
is defined as the conformal mapping of a strip of a two-dimensional (2-D) graphitic lattice 
onto the surface of a cylinder. Each nanotube is uniquely characterized by a 2-D lattice 
vector L = uiSli + n2a2 = [721,^2], where ai,a2 are 2-D primitive lattice vectors. The set 
of integers [ni,n2] determines the geometric properties of the nanotube. The radius of the 



tube is i? = ^ = ^^^^Jnl + 712 + nin2, where do is the C-C bond length. There are two 
symmetry operations: C^ and S{a, h). Cjq is a A^-fold rotation along the axis, where A^ is 
the largest common denominator of Ui and n2- The screw translation S{h,a) represents a 
rotation of a about the axis followed by a translation of h along the axis. The parameters 
h and a are determined from a 2-D lattice vector P = piai +p2Si2, where pi,p2 are integers 
which satisfy the condition p2^i —Pin2 = N ^. In the London approximation, the hopping 
between site i and site j is modified by a phase factor due to the presence of a magnetic field, 
Vij = Voexp (i^ // A(r) -dr). Here Vq is the nearest-neighbor hopping amplitude, A(r) is 
the vector potential associated with the magnetic field, and 0o = ~ is the flux quantum. 

For the special case of a uniform magnetic field H parallel to the tube axis z, both Cat 
and S{a, h) remain symmetry operations of the tube. In this case the hamiltonian can be 
solved analytically and we obtain the following band structure. 



en{K) = ±Vo [3 + 2 cos((5i) + 2 cos(52) + 2 cos((5i + ^2)] 

n = 0, 1, ■ ■ ■ , A^ - 1 , 
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Here k is the pseudo-momentuni associate with the screw translation S{a,h). For H = 
we recover the result of White et al. 0. 

For an idea nanotube each vr orbital contributes one electron, the Fermi energy is at the 
center of the band, ep = 0. From Eq.(l) one finds that, in addition to the dependence on 
the helicity and the radius of the nanotube, the band gap A = 2min{e„(fi;)} varies strongly 
with the magnetic field. Our calculations show that A is a simple period function of the 
fiux, = ttR^H, with a period of (po, 
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An^ 
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30 



for rii — 77-2 = 3g, and 



A{H) 



for rii — 77-2 = 3g ± 1, where 
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is a characteristic energy scale associated with each nanotube. 

From above equations one can draw several conclusions. 1) In the absence of the magnetic 
field there are two types of nanotubes: those with rii — n2 = 3q are metallic and those with 
rii —n2 = 3g±l are semiconducting. 2) In the presence of a magnetic field along the tube axis 
a gap opens up for the metallic tubes. In contrast, the gap decreases for the semiconducting 
tubes, reaching zero at one-third of the fiux quantum. Therefore a field-induced metal- 
insulator transition is expected for all nanotubes. 3) The gap scales inversely with the 
tube radius. Its dependence on the magnetic field is controlled by the characteristic field 
H* = 0o/(vri?^), which can be small for typical nanotubes. For example, taking Vq ~ 2.6eV, 



do = 1.43 A, one obtains A ~ SOmeV for a nanotube of radius R ~ 12nni. For such a gap, 
a magnetic field of 3 T will reduce the gap to zero. Such a field is well within the reach 
of most experiments. 4) Because of the sensitive dependence of the band structure on the 
magnetic field, a large magnetoresistance is expected for all carbon nanotube. 

For a uniform magnetic field making an arbitrary angle 6 with the tube axis, neither Cn 
nor S{a, h) are symmetry operations of the nanotube. But a translation along the tube axis 
with the lattice vector T = ^"^^^"^ ai — "•^J;^"^ a2 is still a good symmetry operation. For this 
case, one can obtain the band structure numerically. It is found that the gap is no longer 
a simple periodic function of the flux. However, our calculations show that for each family 
of nanotubes the reduced gap A{H,6)/Aq is a universal function, when the magnetic field 
is scaled by H*. This universal function is shown in Fig.l for several field directions and 
for both types of nanotubes. The universal scaling relation enables one to estimate the gap 
of any nanotube in an arbitrary magnetic field. From Fig.l one observes that the effect of 
field on the band gap is reduced when the field direction deviates from the tube axis. This 
does not imply that physical properties such as the susceptibility are less affected, as we will 
show. 

The strong field dependence of the band structure suggests a large magnetic susceptibility 
X- There are two contributions to x'- the Pauli term xp (due to the electron spin) and the 
orbital term Xorb (due to the change in the band energy). Because of the low density of 



states near the Fermi energy xp is negligible [O], X = Xp + Xorb ~ Xorb- At T = 0, x can 



be calculated from the second derivative of the total energy, but at a finite temperature the 
free energy should be used |p!3[ , 

d^F{H,T) 



x = - 



F{H,T) = -knTY^ln 



(5) 

H=0 

( en{K,H) - fiY 
1 + exp — 



The summation is over the complete band. It is important to include the effect of the 
magnetic field on the complete band structure because x is calculated from the total energy 
which includes contributions far away from the Fermi level. Thus, a perturbative calculation 



near the Fermi energy gives an inaccurate result in general 

In the case of ideal carbon nanotubes at T = 0, the Fermi energy is at the center of 
the band ep = 0. It is found that for H || z the metallic nanotubes {rii — 722 = 3g) are 
paramagnetic while the semiconducting nanotubes {rii — ^2 = 3g ± 1) are diamagnetic. 
In contrast, for H || z all nanotubes are diamagnetic (Fig. 2 and Fig.3). The magnitude 
of the susceptibility increases linearly with the size of the nanotube regardless whether 
it is diamagnetic or paramagnetic, |x| ~ -R- In general x is a sensitive function of 9. 
Empirically it is found that x(^) = o + bcos{26) is a good approximation (Fig.2). This 
unusual dependence on the field direction suggests that magnetic poling is a possible method 
of aligning nanotubes. 

An important result found is that x is sensitive to the position of the Fermi energy. A 
small change in the carrier density (hence the Fermi energy) from the half filling, which is 
likely in real materials, can lead to a large change in x- The variation of x on ep depends on 
the characteristic energy scale Aq. It is found that x/R as a function of e//Ao is universal for 
each family of nanotubes. Fig.3 shows examples for both H || z and H ± z. One finds that 
in the case of a perpendicular field, x remains diamagnetic in the vicinity of the half filling, 
while for a parallel field a small deviation from the half filling changes the susceptibility 
dramatically. (Notice that for the semiconducting nanotubes, the Fermi energy jumps from 
ei? = to {epl > O.SAq for a small change in the carrier density. For a typical nanotube of 
radius R ~ 12nm, the carrier density at ep = ±Ao is around 1.0 ± 2.5 x 10^^.) 

Finally, in Fig.4 we show the temperature dependence of x- In all cases the magnitude 
of X decreases with increasing T. Again, universal scaling is obtained if the temperature 
is scaled by T* = Ao/ks- For a typical nanotube of radius R ~ 12nm, A ~ 30meV, 
T* ~ 330K. 

Our calculation provides a qualitative explanation for the unusual transport and mag- 
netic properties observed in recent experiments. The large magnetoresistance observed in 
nanotube bundles ||^ is likely related to the field dependence of the band structure near the 
the Fermi level. If so, our calculations predict that the magnetoresistance should depend 

6 



on the field direction and be sensitive to the carrier density. The magnetic susceptibihty 
was measured by two groups, one group for randomly oriented carbon nanotubes 0] and 
the other aligned nanotube bundles [^]. At low temperature unusually large diamagnetic 
susceptibility is found, x ~ —200 to —300 xlO~^cgs/mole for the magnetic field both per- 
pendicular and parallel to the tube axis. From this data we can estimate that the typical 
radius of nanotubes in those samples is around i? ~ 7 to 15nm, a value close to that re- 
ported by experiments. In addition, both experiments found that |x| decreases substantially 
from low temperature to room temperature [Q, suggesting a characteristic temperature 
around SOOii'. From our calculation, this observation corresponds to a typical nanotube 
size of i? ~ 12nm, consistent with the above estimate. One discrepancy remains a puzzle. 
Experimentally, x is niore diamagnetic when the field is parallel to the tube axis than when 
it is perpendicular. Our calculation suggests the reverse in most cases (Fig.3). Though it 
is possible to chose a certain ep such that our results agree with experiments, we feel that 
existing experiments are not well enough controlled for such a quantitative comparison. For 
example, the issue of uniformity of nanotubes and the fact that most of them contain multi- 
ple shells should be addressed. Clearly, further experiments, such as the doping dependence, 
are needed for a quantitative test of the present calculations. 

In conclusion, we have shown that novel properties are to be expected for carbon nan- 
otubes in a magnetic field. A field-induced metal-insulator transition is predicted for all pure 
nanotubes. In general, a large magnetoresistance is expected due to the sensitive dependence 
of the electronic structure on the external magnetic field. The weak-field magnetic suscep- 
tibility is large and increases with the nanotube radius. The susceptibility can be either 
diamagnetic or paramagnetic depending sensitively on the helicity of the nanotube, the field 
direction, and the position of the Fermi energy. A characteristic energy Aq exists for each 
nanotube. Universal scaling is found for x as a function of the reduced Fermi energy ep/^o 
and the reduced temperature /cbT/Aq. Both the magnitude and the temperature dependence 
of the susceptibility are in agreement with recent experiments. Using the current experimen- 
tal data, it is estimated that typical nanotube radii are around i? ~ 7 — 15nm, consistent 



with direct observations. Our results indicate that careful measurements of the magnetic 
susceptibility provides an efficient probe to characterize nanotubes. The directional depen- 
dence of the magnetic response suggests the magnetic poling as a possible method of aligning 
nanotubes. The novel and unusual magnetic properties of nanotubes may have promising 
applications in areas such as the magneto-electronics and the magnetic field detectors. 
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FIG. 1. The reduced gap A/Aq as a function of the magnetic field H for different field di- 
rections. Solid lines: semiconducting tubes (ni — n2 = 3g it 1), ^ = 7r/2, vr/S, 7r/6, (top down). 
Dashed lines: metallic tubes (ni — n2 = 3g), 6 = 7r/2, vr/S, 7r/6, (bottom up). 6 is the angle 
between the field direction and the tube axis, Aq = VqcIq/R. For all figures presented in this paper 
we use Vq = 2.6eV, do = 1.43A. 
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FIG. 2. The scaled susceptibility x/R ^s a functions of the angle 9 between the magnetic field 
direction and the tube axis, ei? = 0, T = 0. Plus: rii — n2 = 3. Diamonds: rii — n2 = 3^ + 1. Lines 
are fits to the functional form a + bcos{26). The case of rii — n2 = 3g — 1 is very similar to that of 
ni — n2 = <iq + I- x is in units of 10~®cgs/mole, R is in units of A. 
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FIG. 3. Universal scaling of x/i? as a function of the reduced Fermi energy epj IS.^. T = 0. Solid 
line: H || z, ni — n2 = 3g. Dotted line: H || z, ni — n2 = 3^ + 1. Dashed line: H _L z, ni — n2 = ?>q. 
Dot-Dashed line: H _L z, rii — n2 = 3(7 + 1. The case of rii — n2 = 3(7 — 1 is very similar to that of 
^1 — ^2 = 3g + 1. xjR is in units of 10~^cgs/mole/A. 
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FIG. 4. The universal dependence of the scaled susceptibility x/^ on the reduced temperature 
ksT/ Aq. The Fermi energy is at ep = 0.2Ao. Squares: H || z,ni — 77-2 = 3q. Diamonds: 
H II z, ni — n2 = 3q — 1. Plus: H ± z, ni — 722 = 3q. Crosses: H ± z, rii — 71,2 = 8(7 — 1. The case of 
^1 — ^2 = 3g + 1 is very similar to that of rii — n2 = 3g — 1. x/R is in units of 10~^cgs/mole/A. 
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